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A detailed account of a recently developed method [Marom et al, Phys. Rev. Lett. 105, 046801 
(2010)] to quantify the registry mismatch in layered materials is presented. The registry index, 
which was originally defined for planar hexagonal boron-nitride, is extended to treat graphitic sys- 
tems and generalized to describe multi-layered nanotubes. It is shown that using simple geometric 
considerations it is possible to capture the complex physical features of interlayer sliding in layered 
materials. The intuitive nature of the presented model and the efficiency of the related computations 
suggest that the method can be used as a powerful characterization tool for interlayer interactions 
in complex layered systems. 



INTRODUCTION 

Single layers of atomically thin molecular structures 
have attracted vast attention in the past few years. Most 
commonly, they appear in quasi-one dimensional [H-[lo| 
and quasi- two dimensional lH-[l8| forms. The phys- 



ical properties of these unique structures depend on 
their specific geometry and chemical c omp osition be- 
ing organicju 0, [l0| inorganic, [M, E E3, Efil or 
mixed. [19l430| When stacked together to form multi- 
layered structures, the physical properties of the in- 
dividual layers m ay b e considerably altered via inter- 
layer interactions. [3ll445j Due to the different nature of 
the intra- and inter-layer interactions, the resulting lay- 
ered systems often present highly anisotropic properties. 
Within each layer covalent bonding usually results in rel- 
atively high strength and stiffness (il-Hol l46| along the 
direction of the layer and, in some cases, in efficient elec- 
tronic (il-Hll. l47j and heat 48Mt3oj transport. In contrast, 
between the layers long-range dispersion and electrostatic 
interactions produce weaker and more flexible mechanical 
and result in reduced transport capa- 



61 



properties 
bilities. 

An important factor governing the physical properties 
of multi-layered materials is the registry matching be- 
tween the layers. Depending on the nature of the in- 
terlayer interactions, different layered materials present 
different optimal stacking modes. As an example, the po- 
lar nature of the interlayer B-N covalent bond in hexago- 
nal boron-nitride (/i-BN) results in considerable Coulomb 
interactions between atomic sites belonging to different 
layers. These, in turn, dictate an optimal A A! stacking 
mode, where a boron(nitrogen) atom in one layer resides 
above a nitrogen(boron) atom in its adjacent layers (see 
Fig. [T]). Graphite is iso-electronic to fo-BN and has the 
same hexagonal structure within each layer. It is there- 
fore tempting to assume that both systems present a sim- 
ilar optimal stacking mode. Nevertheless, due to the lack 
of bond polarization in the homo-nuclear intralayer co- 
valent bonds of graphene the dominant interlayer inter- 
actions in graphite and in few-layered graphene (FLG) 
are dispersion forces. These, in turn, lead to an optimal 
AB stacking mode, where a carbon atom in one graphene 
layer resides on top of a hollow site of the corresponding 



adjacent layers. The picture becomes even more compli- 
cated when considering mult i- walled nanotubes. In such 
systems, apart from the specific chemical composition of 
the nanotube, which dictates the nature of the interlayer 
interactions, curvature differences and different rolling 
orientations (chiral angles) of adjacent layers result in 
complex registry matching and mismatching patterns of- 
ten regarded as Moire fringes, [lif 

It is therefore clear that registry matching plays an 
important role in dictating the electronic properties of 
layered materials and in the interlayer sliding physics of 
such systems. 35, 36|, |45, Nevertheless, previous 

studies concerning such effects have regarded the reg- 
istry matching in qualitative terms labeling it as "bad" , 
"good" , or "optimal" according to the relative energetic 
stability of the different stacking modes as calculated by 
either force fields [62|, [78| or via electronic structure the- 
ories, m, q 

Recently, a quantitative measure of the registry match- 
ing in planar fo-BN was proposed. [45] Based on intuitive 
geometrical considerations and common knowledge re- 
garding the nature of the interlayer interactions in /i-BN, 
a simple model was derived, which predicts the relative 
stability of different stacking modes. It was shown that 
the main features of the interlayer sliding physics in this 
material can be captured by this simplified model thus 
shedding light on the main factors that govern this com- 
plex process. It is the purpose of this paper to give a 
detailed account of this registry index model, extending 
it to other planar layered materials such as graphite, and 
generalizing its applicability to multi-walled nanotube 
structures. 



REGISTRY INDEX IN /i-BN 

In order to define a quantitative measure of the registry 
mismatch in planar /i-BN it is important to understand 
the nature of the interlayer interactions in this system. 
Three important forces should be taken into account: 

1. Dispersion interactions: Van der Waals-London 
forces play a central role in the physics of molecular 
stacking. While being weaker than the intralayer 
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covalent bonding, these induced dipole-dipole in- 
teractions are responsible for anchoring the differ- 
ent layers of a multi-layered material at the appro- 
priate interlayer distance. [IBj Nevertheless, it was 
recently shown that dispersion forces are relatively 
insensitive to the specific interlayer arrangement of 
different stacking modes in /i-BN and thus have lit- 
tle effect on the interlayer sliding potential. [45| 

2. Electrostatic interactions between ionic cores: The 
dominant factor governing the interlayer sliding po- 
tential of /i-BN are electrostatic attractions and re- 
pulsions between the partially charged atomic cen- 
ters on adjacent layers. Due to the difference in 
electronegativity of the two atoms, the boron bears 
a partial positive charge whereas the nitrogen has 



a partial negative charge. |35|, |45|, |80| Based on 
these observations and on basic electrostatic con- 
siderations, one may deduce that optimal registry 
is achieved at the A A' mode which maximizes the 
interlayer N-B Coulomb attractions and minimizes 
the corresponding B-B and N-N repulsions (see 
upper left panel of Fig. [I]). Similarly, the worst 
stacking mode is the AA mode, where the /i-BN 
sheets are completely eclipsed and Coulomb repul- 
sions between atomic centers are maximal (upper 
right panel of the figure). 

3. Charge densities overlap: Another factor that may 
influence the relative stability of different stack- 
ing modes is the electrostatic and Pauli repul- 
sion due to partial overlap of the electron den- 
sities surrounding the boron and nitrogen atomic 
centers. [lH, HH, [80|, EH The electron cloud associ- 
ated with the boron atom is smaller than that of 
the nitrogen atom, one may therefore expect that 
the AB\ stacking mode with eclipsed boron atoms 
(lower left panel of Fig. [I]), will be more energeti- 
cally favorable than the corresponding AB2 mode 
with eclipsed nitrogen atoms (lower right panel of 
the figure). Interestingly, though the interlayer 
electron densities overlap is small [ioj it seems to 
have a considerable effect on the relative stability 
of the AB X and AB 2 stacking modes. [H, EH, [Hj 

Having identified the main interactions involved in h- 
BN interlayer coupling we may now turn to define a quan- 
titative measure of the registry matching in this system. 
Similar to the total energy of the system, we are inter- 
ested in a simple numerical index which will obtain a 
minimum value for the optimal AA' stacking mode and 
a maximum value for the worst AA mode. To this end, 
we ascribe to each atom in the unit cell a circle centered 
around its position (see Fig. [2]). Focusing on two adjacent 
layers, we see that the projection of a circle assigned to 
a specific atom located on one of the layers may overlap 
with circles associated with atoms of the same and/or 
opposite types on the other layer. We mark by Sij the 
overlaps between two such circles, one associated with 
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FIG. 1: High symmetry stacking modes of /i-BN. Upper left 
panel: the optimal AA! stacking mode. Upper right panel: 
the fully eclipsed AA stacking mode. Lower left panel: the 
meta-stable AB\ stacking mode. Lower right panel: the high 
energy AB2 stacking mode. Blue(yellow) circles in the on-line 
version represent boron (nitrogen) atoms. Dashed gray lines 
with large circles represents a lower /i-BN layer and solid black 
lines with small circles represents an upper layer. 



an i atom on one layer and the other with a j atom on 
the second layer, i and j being either B or N. It is now 
clear that the sum Snn + Sbb — Snb complies with our 
requirement of obtaining a minimum(maximum) value 
at the AA' {AA) stacking mode, where Snb is maxi- 
mal(minimal) and Sbb and Snn are minimal (maximal). 
By normalizing this sum to be limited to the range [0,1] 
we obtain the registry index {RI) for /i-BN: 
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Here, S^ A and S^ A are the corresponding overlaps at 
the A A' and A A stacking modes, respectively. 

This definition was obtained based on the knowledge 
that in /i-BN perfect (worst) registry is achieved at the 
AA' {AA) stacking mode due to electrostatic interactions 
between partially charged atomic centers. As stated 
above, the effect of interlayer overlap of charge densities 
influences the relative stability of the AB\ and AB2 con- 
figurations. In order to take this into account we assign 
different radii to circles associated with boron (re) and 
nitrogen (r/v) atoms. This may be viewed as a simplified 
representation of the difference in atomic radii between 
the partially negatively charged nitrogen and positively 
charged boron atoms. For simplicity, we fix the circle 
associated with the nitrogen atom to half the B-N bond 
length, and use the ratio ^ as a free parameter. By 
choosing — < 1, RI obtains a lower value for the AB\ 
stacking mode with respect to the AB2 mode thus repro- 
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FIG. 2: Registry index definition of the overlap area between 
circles assigned to atomic positions in the upper layer (trans- 
parent circles) and their lower layer counterparts (opaque cir- 
cles). The circles representing the atomic centers in Fig. [1] 
were omitted for clarity. The color code is the same as in 

Fig.m 



during the physical requirement. 

It is now possible to plot the RI at different stacking 
modes and compare the resulting surface obtained from 
simple geometric considerations with the sliding energy 
surface obtained from advance electronic structure calcu- 
lations. Such a comparison was recently presented, [45| 
showing a remarkable agreement between density func- 
tional theory (DFT) results obtain via the PBE den- 
sity functional approximation [83| augmented with the 
Tkatchenko- Schemer Van der Waals correction 0, EHJ 
and the RI model with ^ = 0.3. This exemplifies 
the fact that the /i-BN sliding process is governed by 
registry mismatch via electrostatic interactions and val- 
idates our assumption regarding the choice of different 
boron and nitrogen circle radii within the RI model. We 
therefore conclude that the RI, as defined above, can be 
used to characterize the different stacking modes of h- 
BN. The question arises whether this simple geometric 
model can be extended to treat other layered materi- 
als such as graphite, and generalized to more complex 
structures such as nanotubes. Should the answer to this 
question be positive, one would be able to gain valuable 
physical intuition regarding such layered materials and 
characterize their relative interlayer configurations at a 
fraction of the computational cost of current electronic 
structure methods. 

A clue to the answer to this question can be found 
in related recent studies using geometric considerations 
for the description of halogen atoms and rare gases ad- 
sorbed on (111) metal surfaces. 8(3, 8^ In what follows 
we show how the current RI model can be extended and 
generalized as suggested above. 



REGISTRY INDEX IN GRAPHITIC MATERIALS 

We start by showing that the RI model is not lim- 
ited to the case of /i-BN and can be extended to treat 
other planar layered materials. We consider the case of 
graphite or FLG. As stated above, because of the homo- 
nuclear nature of the bonds in these systems no charge 
polarization occurs. Therefore, the main factors govern- 
ing the sliding physics are dispersion forces and overlap 
of charge densities. Due to the lack of electrostatic forces, 
the A A stacking mode, which is equivalent to the stable 
A A' mode in /i-BN, is found to be a maximum on the in- 
terlayer potential energy surface. Furthermore, the AB 
configuration, which minimizes charge densities overlap, 
is the optimal stacking mode of graphite. Interestingly, 
the interlayer distance, which is mostly influenced by dis- 
persion interactions, [45|, l88l - [93| is found to be very 
similar to that of h-BN (3.35 A vs. 3.33 A in ft-BN). 

Based on these observations one can now extend the 
registry index definition to treat graphitic materials. The 
graphitic RI should obtain a minimum value for the op- 
timal AB stacking mode of graphene and a maximum 
value for the worst AA mode of this system. As for the 
case of /i-BN, we ascribe to each atom in the unit cell a 
circle of radius rc = 0.5d cc centered around its position, 
where d cc is the intralayer carbon-carbon covalent bond 
length. The overlap between two such circles, one associ- 
ated with a carbon atom on one layer and the other with 
a carbon atom on the second layer, is then marked by 
Sec- Naturally, if the RI is chosen to be proportional 
to Sec it will comply with the requirement of obtaining 
a minimum(maximum) value at the AB(AA) stacking 
mode. By normalizing RI to be limited to the range 
[0,1] we obtain the following definition: 



_ Scc-S$S 
^graphitic gAA _ gAB ' 



(2) 



Here, Sqq and Sqq are the overlaps at the A A and AB 
stacking modes, respectively. 




FIG. 3: Registry index surface of double layered graphene. 
In Fig. [3] the registry index surface is presented as a 



function of relative interlayer sliding parallel to the basal 
planes of a bilayer graphene system. Comparing to re- 
cent molecular dynamics [78| a nd dispersion-augmented 
tight-binding calculations, [79| it is found that the RI 
landscape fully captures all the important features of 
the interlayer sliding physics of graphene. Furthermore, 
our results are consistent with recent experimental and 
theoretical investigations showing an orientation depen- 
dent sliding resistance in graphitic systems. @,[95| This 
proves that the RI concept is of general nature and can 
be readily extended to characterize the interlayer inter- 
actions in a variety of layered materials. We now turn 
to describe how this model can be further generalized to 
treat more complex structures. 



REGISTRY INDEX IN MULTI- WALLED 
NANOTUBES 





FIG. 4: Schematic representation of the procedure to calcu- 
late the RI of DWNT systems. First, the two layers are un- 
rolled. Next, the narrower sheet (inner layer) is stretched to 
match the width of the wider sheet. Finally, circles are placed 
around the atomic positions and the RI surface is calculated 
for different relative positions of the unrolled layers. These, 
in turn, are equivalent to relative telescoping and rotation of 
the tubular system. 



Planar layered materials usually have a compact unit 
cell, which can be readily treated using standard elec- 
tronic structure methods within periodic boundary con- 
ditions calculations. On the other hand, despite their 
reduced dimensionality, even achiral nanotubes present 
relatively large unit cells. This is especially true in the 
case of multi-walled nanotubes (MWNTs) where often 
one finds that, apart from the smallest bilayer systems, 
they are beyond the reach of state-of-the-art electronic 
structure methods. It is therefore desirable to generalize 
the RI defined above to treat tubular structures. Once 
such a generalization is established, it can be used as an 
efficient and reliable characterization tool for the relative 
stability of different inter-tube configurations. 

Since nearest neighboring layers interactions are the 
most important factors governing the relative arrange- 
ment of the different layers within a MWNT, we gen- 
eralize the RI to the case of double-walled nanotubes 
(DWNT). This allows the investigation of the isolated 
layer-layer interactions, which are at the basis of the 
multi-layered system behavior. 

The procedure to calculate the RI in DWNTs is 
schematically presented in Fig. HI We start by cutting 
the two layers along a given line parallel to the axis of the 
tube. Next the layers are unrolled to form planar sheets 
of different width. Then, the narrower sheet (unrolled 
inner tube) is stretched to match the width of the wider 
sheet (unrolled outer tube), thus taking into account the 
effect of curvature on the registry mismatch between the 
two layers. Finally, circles are placed around the atomic 
centers of the two layers and the RI is calculated using 
Eq. [D (or Eq. [2]) for different interlayer shifts parallel to 
the basal planes of the two layers. The resulting RI sur- 
face corresponds to relative telescoping and rotation of 
the two tubes within the DWNT. 

Similar to the case of fo-BN, we can now compare 
the RI surfaces to the results obtained by DFT calcula- 
tions. To this end, we perform a set of DFT calculations 
with the Gaussian suite of programs. [§6|, [97| The lo- 



cal spin density (LSDA), PBE, [H and HSE06 H ^ 
exchange-correlation functional approximations are used 
together with the double-^ polarized 6-31G** Gaussian 
basis set jlOOj . As discussed above, Dispersion interac- 
tions play a major role in anchoring the layers of fo-BN at 
the appropriate interlayer distance. Nevertheless, in the 
case of DWNTs, the interlayer distance is fixed by the 
differences of curvature between the two tubes which are 
set by the tubes indices. Hence, the effects of dispersion 
interactions on the interlayer sliding energy, which have 
been shown to be of minor importance in fo-BN, [45| are 
neglected in the present work. 

Results for three representative double-walled 
boron-nitride nanotubes (DWBNNT) are presented: 
(5,5)@(10,10), (6,6)0(11,11), and (6, 0)@(14, 0), where 
the notation (ni, rai)@(ri2, 7712) stands for an inner 
(ni,rai) tube placed within an outer (n2,ra2) tube. 
Unlike the case of planar fo-BN, DWNTs present a wide 
range of possible structures. The two tubes may differ 
in chiralities, a factor that may considerably alter their 
registry matching and result in orders-of-magnitude 
differences in their sliding energy corrugation. Interest- 
ingly, even for achiral tubes of the same kind (armchair 
or zigzag) two types of systems can be constructed: the 
first (type-I) resulting from rolling two h-BN sheets in 
the AA! stacking mode, and the other (type-II) resulting 
from rolling two A A stacked fo-BN layers. It should be 
noted that one such achiral DWNT may be obtained 
from its counterpart by switching the identities of the 
boron and nitrogen atoms in one of the layers. As shown 
below, once the chiralities and types of the two tubes are 
set, changing diameters of the tubes, even while fixing 
the inter-tube distance, has remarkable impact on the 
registry matching between the layers and their sliding 
energy surface corrugation. 

In Fig. [5] the results for the type-I (5, 5)@(10, 10) sys- 
tem are presented. The tubes are formed by rolling two 
A A! stacked layers, while fixing the B-N distance to be 
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FIG. 5: Rotation-telescoping energy landscape of the type-I 
(5,5)@(10,10) armchair DWBNNT. Upper left panel: Relative 
total energies of different inter-tube configurations calculated 
using DFT at the PBE/6-31G** level of theory. Lower right 
panel: Registry index surface calculated using the procedure 
described in the text. Axial and side views of the system are 
shown to emphasize the effects of curvature on the registry 
mismatch between the two tubes. 



~1.44 A. No geometry optimizations are performed. The 
resulting distance between the tubes is ~3.44 A which is 
similar t o the equilibrium interlayer distance of fo-BN of 
3.33 A. [lOl| In the upper left panel of the figure re- 
sults obtained at the PBE/6-31G** level of theory are 
presented. Similar results have been obtained using the 
LSDA and HSE06 functionals (not shown). The inter- 
layer potential energy is found to be much more sensitive 
to relative rotations of the two armchair tubes than to 
telescoping parallel to the tubes axis. The corrugation 
energy, which is defined as the maximal amplitude of 
energy changes between different interlayer relative po- 
sitions, is found to be ~ 0.02 eV/unit-cell. The corre- 
sponding RI surface presented in the lower right panel 
of the figure reproduces all of these effects while cap- 
turing even the fine details of the sliding energy surface 
landscape. 

Fig. [6] presents similar results, obtained at the 
LSDA/6-31G** level of theory, for the type-I 
(6,0)@(14,0) zigzag DWNT system. As can be 
seen from the upper left panel, the sliding energy surface 
is very similar to that obtained for the (5,5)@(10,10) 
armchair system. Two important differences are ap- 
parent: (i) the role of the axes is interchanged (ii) 
the corrugation energy for the zigzag system is found 
to be an order of magnitude larger than that of the 
armchair system. The reason for the latter difference is 
the smaller interlayer distance of 3.18 A in the zigzag 
DWNT system. Nevertheless, as in the case of the 
armchair DWNT, the RI landscape fully captures all 
the details of the sliding energy surface obtained via 
DFT calculations. 

To exemplify the complexity of the sliding physics of 
DWBNNT systems, [El the (6,6)@(11, 11) is considered 
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FIG. 6: Rotation-telescoping energy landscape of the type-I 
(6,0)@(14,0) zigzag DWBNNT. Upper panel: Relative total 
energies of different inter-tube configurations calculated using 
DFT at the LSD A/6-31G** level of theory. Lower panel: Reg- 
istry index surface calculated using the procedure described 
in the text. 



as well. This armchair DWNT has the same interlayer 
distance as the (5, 5)@(10, 10) system considered above. 
One may naively expect that the sliding energy land- 
scape of the two systems, which have the same chiral- 
ity, type, and interlayer distance, would be similar. As 
can be seen in the left panel of Fig. [7J this is not the 
case. The Corrugation energy of the (6, 6)@(11, 11) is 
found to be an order of magnitude smaller than that of 
its (5, 5)@(10, 10) counterpart. In fact, the energy dif- 
ferences between relative tube positions are calculated 
to be smaller than 0.002 eV/unti-cell, which is beyond 
the accuracy of our DFT calculations. Accordingly, the 
agreement between the RI landscape (right panel of the 
figure) and the DFT results is not as good as those ob- 
tained for the (5,5)@(10, 10) and (6,0)@(14,0) systems. 
Consistent with the reduction in the corrugation energy, 
the magnitude of the RI variations is reduces by more 
than an order of magnitude as well. It should be stated 
that the RI remains a valid quantity to describe the reg- 
istry matching in this system. Furthermore, the ampli- 
tude of the RI variations may serve as an indication to 
the ability of DFT calculations to accurately describe 
the interlayer sliding landscape. It remains to be shown 
whether in such cases of extremely small corrugation en- 
ergy the sliding physics is still dominated by the reg- 
istry mismatch. An answer to this question can be given 
only with more accurate electronic structure calculations 
including the detailed effects of dispersion interactions, 




which may have an important contribution in these situ- 
ations. 
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FIG. 8: Comparison of the recurrence frequency of on- 
top B-N stacking arrangements in type-I (5,5) @( 10, 10) and 
(6,6)0(11,11) DWBNNTs. 



tity of its layers. Nevertheless, the extension of the RI 
model to treat tubular structures proves to be a reliable 
tool for the quantification of the registry matching be- 
tween the layers. Therefore, it can be used to character- 
ize the interlayer potential and identify optimal interlayer 
configurations of very large multi- walled nanotubes that 
are beyond the reach of current DFT calculations. 



FIG. 7: Rotation-telescoping energy landscape of the type- 
I (6,6)0(11,11) DWBNNT. Upper panel: Relative total en- 
ergies of different inter-tube configurations calculated using 
DFT at the LSDA/6-31G** level of theory. Lower panel: Reg- 
istry index surface calculated using the procedure described 
in the text. 

In order to better understand these differences in 
the corrugation energy and RI variations between the 
(5,5)0(10,10) and (6,6)0(11,11) DWBNNT we take a 
closer look at their symmetry characteristics. We choose 
an inter-tube arrangement which has a perfect on-top 
stacking between a given boron atom on one wall and a 
nitride atom on the other wall (see white lines in Fig.[8j). 
It is now possible to define a recurrence frequency as 
the number of times such an on-top stacking appears 
along the tube circumference at that given geometry. 
Since a (n, n) boron-nitride nanotube has a n-fold ro- 
tational symmetry around the axis of the tube, the re- 
currence frequency of a (ni, ni)0(ri2, 712) system is given 
by gcd(n\,ri2), where gcd stands for the greatest com- 
mon divisor. For the two armchair DWNTs considered 
above, the recurrence frequencies are gcd(5, 10) = 5 and 
0cd(6,ll) = 1 for the (5,5)0(10,10) and (6,6)0(11,11) 
systems, respectively. Naturally, as the recurrence fre- 
quency grows, the RI (and the total energy) of the 
on-top configuration decreases, and the corrugation of 
the sliding energy surface increases, explaining why the 
(5, 5)0(10, 10) system presents a considerably higher cor- 
rugation energy than the (6,6)0(11, 11) DWBNNT. 

Similar considerations can be used to characterize 
DWNT of other chiralities showing the sensitivity of the 
sliding energy surface of a DWBNNT to the specific iden- 



SUMMARY 

A new methodology to quantify the registry match- 
ing in layered materials, based on simple geometric con- 
siderations, was presented. The registry index, which 
was originally developed to describe the stacking reg- 
istry in planar /i-BN, was extended to treat graphitic 
materials and generalized to describe multi-walled nan- 
otubes. Even in the challenging case of double- walled 
boron-nitride nanotubes, the RI model was able to cap- 
ture the important physical features of the interlayer slid- 
ing up to fine details. This marks the method as a pow- 
erful characterization tool for interlayer interactions in 
complex layered systems while giving intuitive insights 
regarding the nature of the interlayer couplings. 
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